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Abstract 

We consider a lattice model for three-component systems in which 
the lattice bonds are covered by molecules of type AA, BB, and AB, 
and the only interactions are between the molecular ends of a common 
lattice site. The model is equivalent to the standard Ising model, and 
the coexistence curves for different lattices and/or some specific form 
of the interactions have been previously investigated. We derive the 
spinodal curve of the three-component model on the honeycomb lattice 
based on the mean-field and Bethe-lattice results of the equivalent 
Ising model. The spinodal and the coexistence curves of the ternary 
solution are drawn at different values of the reduced temperature, the 
only parameter of the model. The particular case of a two-component 
system is also illustrated. 
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1 Introduction 



Two almost completely immiscible liquids may increase their mutual sol- 
ubility in the presence of a solvent. By increasing the solvent concentration, 
the two liquid layers become more alike and, at some (temperature depen- 
dent) composition, a single homogeneous phase is formed. In order to de- 
scribe this phenomenon, Wheeler and Widom |T[ considered three types of 
diatomic molecules, AA, BB, and AB, covering the bonds of a regular lattice 
(one molecule per bond) and subject to the condition that only A atoms or 
only B atoms may meet at a given lattice site (infinite repulsion between 
A and B atoms of a common site, no interaction otherwise). The model 
reduces to the standard (spin-1/2, nearest-neighbor interaction) Ising model 
on the same lattice and its ferromagnetic transition corresponds to a phase- 
separation transition for the ternary solution. Although the model catches 
the main aspects of the phase-separation transition, it is not temperature de- 
pendent. The model can be extended to include temperature clS du SI gnificant 
variable by allowing finite interactions between any two atoms of a common 
site. This extended Wheeler- Widom model is still equivalent to the standard 
Ising model but on a decorated lattice in which each site of the original lattice 
is replaced by a cluster of sites. The exact coexistence (binodal) curves of 
this extended Wheeler- Widom model on the honeycomb and Bethe lattices 
were determined in || and |§ , respectively. Further extensions of the model 
have been considered in H ^ ||. 

A solution suddenly quenched from an initial one-phase state in the two- 
phase region reaches its equilibrium through either a nucleation or a spinodal 
decomposition process f7j. The states below the coexistence curve can be 
then classified as metastable (close to the binodal) and unstable (deep in the 
two-phase region), respectively. The spinodal curve is supposed to separate 
these two types of behavior. According to the classical (mean-field) theory 
of first-order phase transitions, the spinodal is defined as the locus for which 
an appropriate susceptibility diverges; for a multi-component solution, it 
follows that the determinant of the diffusion coefficients must be zero and 
this provides an experimental determination of the spinodal curve. Although 
a mean-field-like approach oversimplifies the problem, it is still a convenient 
way to distinguish between the metastable and unstable states. 

The aim of the present paper is to determine the spinodal curve for the 
extended Wheeler- Widom model. The molecules are considered to occupy 
the bonds of a honeycomb lattice (see Fig. [I]) and the lattice of the corre- 
sponding Ising model is called 3-12 (see Fig. 0). As shown in Section 2, by 
applying at first a star-triangle transformation and then a double-decoration, 
the Ising model on the 3-12 lattice becomes the Ising model on the honeycomb 
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lattice, the original lattice of the molecular model. In Section 3 we shortly 
review two approximate treatments of the ferromagnetic Ising model: mean- 
field and Bethe-lattice. The resulting coexistence and spinodal curves for 
the three-component model, the effect of temperature, and the limit of a 
two-component system are all presented in Section 4. In the last section we 
discuss some possible developments of the present work. 



2 Extended Wheeler- Widom model on the 
honeycomb lattice 

The grand-canonical partition function of the Wheeler- Widom model 
with finite interactions on the honeycomb lattice (N sites) is proportional to 
the partition function of the standard Ising model on the 3-12 lattice 

-^3-12 = Yl eX P \ R a i a i + L Yl ° i0 3 + h J2°~i) 

{<Ti} \ (M')CC 3 (i,j)cC 2 i J 

with <7j, (Tj = ±1 (i, j = 1, . . . , N) depending on the spin orientation. The 
first sum in ( |2.1| ) is over all possible configurations, the second one runs over 
all the triangles (C3 graphs) of the 3-12 lattice, the third one is over all the 
bonds (C2 graphs) connecting the triangles, and the last one is over the whole 
lattice (see Fig. ||). The parameters R, L, and h from (|2.1|) are related to 
the parameters of the original model as follows 

R= Ti m (2 £ AB - £aa- ebb), (2.2) 
(Haa + Hbb - Z^ab), (2.3) 



4k B T 
and 

h= ——(^aa- I^bb -2e AA + 2e BB ). (2.4) 
Ak B T 

Here k B is Boltzmann's constant and T is the temperature. We use Eaa, £ab, 
and e BB to denote the interaction energies between the molecular ends of a 
common lattice site. The chemical potentials Haa, Hbb, and \i A b ah tend to 
infinity (no vacant sites are allowed) , but differences of any two of them are 
finite thermodynamic variables. 

The partition function of the standard Ising model on the 3-12 lattice can 
be connected to the partition function of the standard Ising model on the 
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honeycomb lattice 



(2.5) 



via two transformations: star-triangle followed by double-decoration (see 
Fig. |). The result is M 



where 
and 



Zg_ 12 (i2, L, h) = A~ N (Li)B 3N ^ 2 (Li, L, h)Z hon (K, 3h* 
A 4 (Li) = 16cosh(3Li)cosh 3 (L!) 



cosh(2Li) = (e m + l) /2 , i2 > 0, 
are obtained from the star-triangle transformation, and 



(2.6) 

(2.7) 
(2.8) 



B 4 



16e 



-4L 



1 + e 2L cosh(2Lx + 2/i) 1 + e 2L cosh(2Li - 2/i) 



.2L 



AK 



cosh(2Li) + e 2i cosh(2/i) 
1 + e 2L cosh(2L! + 2h)] \l + e 2i cosh(2L! - 2/i) 



and 



[cosh(2Li) + e 2L cosh(2/i)] 2 
4fe 1 + e 2L cosh(2£ a + 2/i) 



(2.9) 
(2.10) 

(2.11) 



1 + e 2L cosh(2Li - 2h) 
follow from the double-decoration transformation. The parameter L\ from 
equations (gl| to O 

is an interaction constant of the Ising model on the 
intermediate lattice (see Fig. |3|); its sign is undetermined and, without loss 
of generality, we may choose L\ > 0. 

Let us denote by Xaa, Xbb, and X^b the mole fractions respectively 
corresponding to the three types of the molecules. The mole fractions obey 
the conservation equation 

Xaa + X BB + X A b = 1 (2-12) 
and they can be related to the quantities of the equivalent Ising model by f2j 

Xaa + X B b — Xab = 03-12 (2-13) 

and 



Xaa — X BB — m 3 _ 12 , 



(2.14) 
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where cr 3 _ 12 =< OiOj >c 2 describes the spin-spin correlation on the C2 graphs 
and m 3 _ l2 =< cr^ > is the magnetization of the 3-12 lattice. Using ( |2.6|) , cr 3 _ 12 
and m 3 _ 12 can be expressed in terms of the similar quantities a and m on the 
honeycomb lattice as 



a 



3-12 



dL 



+ 



Li,h 



'die 

dL 



a 



Li,h 



dh* 
~dL 



m 



Li,h 



(2.15) 



and 



m 



3-12 




Li,L 




a 



'9K 
dh 



m. 



(2.16) 



The Ising model has two parameters, K and h+, while the Wheeler- Widom 
model has three independent combinations of the original parameters: R, L, 
and h. One can choose 1/R as the reduced temperature of the molecular 
model 0; we will use t to denote it and keep it constant throughout the 
calculation. (When eab — * 00, then t — >• and the model reduces to the 
original model introduced by Wheeler and Widom [|TJ; see also Appendix.) 
Because we describe the molecular model by means of the Ising model, all 
the derivatives from the equations ( p. 15 ) and ( 2.16Q must be expressed in 
terms of R, K, and rather than L\, L, and h. The interaction constant L\ 
as a function of R is given by ( |2.8| ) and from ( |2.10| ) and ( |2.11| ) we can find L 
and h as functions of K and at constant L\ (or, equivalently, at constant 
t). Namely, 



tanh(2/i) 



sinh(2Li) sinh(2/ij 
cosh(2Lx) cosh(2/iJ - e - 2K 



X 



and 



,2L 



cosh(2L!) cosh(2^ 



-2K 



= A. 



(2.17) 



(2.18) 



VI -X 2 cosh(2L l )e- 2K - cosh(2^ 

The expressions ( |2.17|) and ( ^.18[ ) are well defined only if |%| < 1 and A > 0, 
which amounts to 

\K\ < L i (2-19) 



and 



- -ln[cosh(2Li -2\h+ 



< K < - In 
2 



cosh(2L! 
cosh(2/i^ 



(2.20) 



the latter inequalities determining the region in the (K, h^-sp&ce which can 
be mapped into (L, h; i?)-space. This region becomes the whole space when 
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R —>■ oo. The derivatives from ( |2.15 ) and (|2.16| ) as functions of K, h*, and 
R are given in Appendix. Provided a and m (both functions of K and h+) 
are known, we have now the complete "dictionary" between our molecular 
system and the Ising model. 

The ferromagnetic transition of the Ising model (K > K c > 0, h+ = 0) 
corresponds to a phase-separation transition for the Wheeler- Widom model 
(L > 0, h — 0). If (T (K) = a(K,0) and the spontaneous magnetization 
m Q (K) = m(K,0) from ( |2 . 1 5| ) and ( |2.16| ) are specified, then one can deter- 
mine the isothermal coexistence curves for our molecular system by elimi- 
nating K between the equations ( 2.13| ) and ( 2.14|) . The coexistence curves 
based on the exact solution of the Ising model on the honeycomb lattice and 
on the three-coordinated Bethe lattice were determined in the and ||, 
respectively. In finding the spinodal curve, we need a s {K) = a[K , h[ s) (K)} 
and m s (K) = m[K,h*(K)] calculated along a certain path h[ s \K) of the 
Ising model. We consider two approximate methods below: mean-field and 
Bethe-lattice. 



3 Ising ferromagnet in the mean-field approx- 
imation and on the Bethe lattice 

In the mean-field approximation, the magnetization of the Ising model 
with the partition function (|2.5|) is the solution of the equation 

m = ta.nh(qKm + 3h+), (3-1) 

where q is the coordination number of the lattice. This equation is invariant 
under the transformation (m, h^) — > (— m, — h^), i.e., is an odd function 
of m. The spin-spin correlation function a can be determined, for example, 
by using its relation to the internal energy per spin u, 

ff = ^(sf + 3ft - m )- (32) 

and substituting the expression for u from || into (|3.2|). As a result, 

a = m 2 . (3.3) 

In zero magnetic field and for K > K c = 1/q the model exhibits a ferro- 
magnetic transition. The spontaneous magnetization is given by ( |3.1| ) with 
= 0: if m is a solution, then — m is also a solution and both these 
values minimize the free energy as function of the magnetization. Once 
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having determined the spontaneous magnetization, the relation (|3.3j ) gives 
o" = m\. At constant K, the "van der Waals loops" of h+ as function of 
m (for |m| < |w |) are interpreted within mean-field-like theories as non- 
equilibrium states: metastable for |m| > |m | and unstable for |m| < |m 



where ±m s give the extremum points of (which, at the same time, are the 
inflexion points of the free energy as a function of the magnetization). By 
changing K, these two extremum points describe the spinodal curve, and the 
mean-field value of m„ is 



s 



/ K 1 
m s = ±y 1 - with K c = -. (3.4) 

The spin-spin correlation function along the spinodal curve is then given by 
a s = m 2 . 

In an improved scheme, with the same results as the Bethe approximation 
10| , we consider the Ising model on the Bethe lattice (see Fig. |j). Defining 



z = e ~ 2K , li = e _6fc *, and ^ = Lix q -\ (3.5) 
the magnetization of the Ising model on the Bethe lattice with the coordina- 



tion number q is [10 



1 + 2zli 1 + [i{ 
while the spin-spin correlation function is [|ll|] 

1 — 2zli 1 + n\ 



m = - ■ - t — (3-6) 



a 



1 + 2zil x + //?' 
where \i x obeys the equation 

fh = f Z + Vi V" 1 
H \1 + zllJ 

From (|3.6|) and ( |3.7|) we can get a as a function of z and m, 



(3.7) 



(3.8) 



1 — m 2 i v 

a = \-2z =. (3.9) 

z + Jl — (1 — z 2 )m 2 

This is the analog of Eq. ( |3.3|) of the mean-field approximation and, for 
z — * 1 (if — > 0, i.e., small couplings), they coincide. In zero magnetic field 
and for z < z c = 1 — 2/q the model undergoes a ferromagnetic transition. 
Both m and a can be determined in principle by solving at first ( |3.8j ) with 



7 



\i = 1 and then substituting the solution fi^z) into ( |3.6| ) and (|3.7|). In the 
particular case of q = 3 one gets || 



1 /1-3A 1/2 



and 



(Tn = 1 + 2Z 



1 - 2z V 1 + 2 
1 1 



(3.10) 



(3-11) 



_l-z 3(1 + z) 3(1 -2z) 

In determining the spinodal curve, we should first derive function of 

m using (|3.6|) and Q3.8|) , and then one readily shows that /i^ has extremum 
points at 

2 



m„ = ±. — ^ with z„ = 1 

\ 1 - z 2 q 

Combining ( |3.12| ) with ( |3.9| ) we easily find a s . 



(3.12) 



4 Phase diagrams 

The ferromagnetic transition of the Ising model occurs in zero magnetic 
field and for a coupling constant K greater than a critical one K c . This tran- 
sition has a correspondent in the molecular model only if K max from ( |2.20| ), 



c 



calculated at = 0, is greater than K c . The condition ln[cosh(2Li)] = 2K t 
then determines the upper consolute temperature t u of the model, i.e., the 
maximum temperature at which a phase separation transition can occur in 
the solution, 

4 

tu = In [2e 2 ^ - 1] ' ( 4-1) 

In the mean-field approximation of the Ising model (q = 3) we get t u ~ 3.762; 
for the Bethe lattice, one has t u = 4/ In 5 ~ 2.485. 

In drawing the isothermal coexistence and spinodal curves, we use the 
triangular diagram: a fixed composition is represented by a point inside an 
equilateral triangle of the unit height, the mole fraction of each component 
being given by the distance from that point to the corresponding side of 
the triangle. For a triangular diagram, the conservation equation ( |2.12| ) is 
automatically fulfilled, and the correspondence with the three components 
of the molecular solution is chosen as shown in Fig. |5]. 

Both the coexistence and the spinodal curves are symmetric under reflec- 
tions around the (Xaa = ^Bs)-axis, a consequence of the symmetry proper- 
ties of the derivatives from (|2.15|) and (|2.16|) when — > —h+ (see Appendix). 
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The analog of the Curie point for a ferromagnet is the plait point; it is given 
by the limit K —>■ K c and it lies always on the (Xaa — ATe^-axis. At 
the plait point, the spinodal curve is tangent to the coexistence curve. In 
the vicinity of the plait point, both the coexistence and the spinodal curves 
can be approximated by parabolas. In the opposite limit, at K = K max , it 
is shown in Appendix that Xab = and, thus, both the coexistence and 
spinodal curves reach the base of the triangular diagram. 

The phase diagrams at some particular values of the reduced tempera- 
tures are presented in Figs. |5] and |], respectively corresponding to the mean- 
field and Bethe-lattice approximations of the equivalent Ising model. In the 
mean-field approximation (Fig. |5|), at t — the spinodal and the coexistence 
curve coincide; they are described by the same parabola y — (1 — 3x 2 )/2. 
This means that no metastable states are predicted at zero temperature, an 
expected result from more general considerations. Technically, this is a con- 
sequence of the relation ( |3.3| ) between a and m, which is independent of the 
coupling constant K: the distinct K- dependencies of m for the spinodal and 
coexistence curves play only a scaling role (the same curve is covered in dif- 
ferent ways). In the Bethe-lattice case (Fig. |6|), the relation (|3.9|) between a 
and m explicitly depends on the coupling constant, and, in general, it should 
be so. Consequently, the three-component system can be in a metastable 
state even at zero temperature, a fact that might reveal an inconsistency of 
the theory. (The implications of this result and a possible way to overcome 
the difficulty within mean-field-like theories of first-order phase transitions 
will be discussed in a forthcoming paper.) By increasing the temperature, 
the curves go down and shrink to a point when t = t u . The limiting case 
Xab = corresponds to a two-component system. In the mean-field ap- 
proximation of the equivalent Ising model, the spinodal and the coexistence 
curves depend on temperature as in Fig. [7]; the case of the Bethe lattice is 
illustrated in Fig. |[ 

5 Discussion 

The Wheeler- Widom model (with finite, two-body interactions) for three- 
component molecular systems is equivalent to the standard (spin-1/2, nearest- 
neighbor interaction) Ising model. Based on this equivalence, we derived the 
spinodal curve for the molecular model on a honeycomb lattice using ap- 
proximate (mean-field and Bethe-lattice) results for the Ising model. Both 
the spinodal and the (previously determined) coexistence curves have been 
drawn at different temperatures. 

The present work can be extended in several ways. One direction is 
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to consider three-body interactions between the molecular ends associated 
with the same site of the honeycomb lattice. The coexistence curves then 
become asymmetric, with the plait point located off the (Xaa = Xg^-axis 
of the triangular diagram ||; a similar asymmetry should be observed for the 
corresponding spinodal curves. Another possible development is to determine 
the spinodals for systems with both upper and lower consolute temperatures. 
Within the present formalism, a lower bound on the transition temperature 
follows from the assumption that the molecular ends may form bonds; the 
coexistence surface for a three-coordinate Bethe lattice has been determined 
in H . A more general situation, with two possible transitions (one with both 
upper and lower consolute temperatures, the second one only with an upper 
consolute temperature below the minimum transition temperature of the first 
one) has been analyzed in ||. It is also possible to generate phase diagrams 
with all the characteristics mentioned above (asymmetric, upper and lower 
consolute temperatures, etc.) within a single ternary solution model of the 
Wheeler- Widom type [Q. 

In the framework of the formalism presented above, the spinodal curve 
of a ternary solution is determined starting from a simple microscopic model 
reducible to two components. Together with analysis of experimental data, 
it could contribute to a better understanding of interesting phenomena such 
as diffusion and nucleation in multi-component systems. Another impor- 
tant aspect should also be mentioned: since a Wheeler- Widom-type model 
is equivalent to the standard Ising model, it perhaps provides the only way 
to get exact results for systems with three components. 
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Appendix 

Let us denote the derivatives that appear in ( 2.15|) and ( 2.16|) by 



Cll 



Cl2 



'd\nB" 
~dL~ \ 

'dK\ 



°2l = 7T 



L lt h 



1 (dlnB 
2 



i 



dL 



Li,h 



Cl3 

Using the notations 



C22 



C23 



2 \ 9/i 



(A.l) 



'dh* 

dh 



Li,L 



-2K 



p = ( e 4R + l) /2 and z = e 

the quantities x an d ^ respectively defined in (|2.17f) and ( 2.18 ) become 
y/p 2 — 1 sinh(2/i + ) pcosh(2/i + ) — £ 



(A.2) 



X 



pcosh(2/iJ — z 



and A 



p,2 — cosh(2/^ 



Introducing 



W = (l+pA) 2 -(p 2 -l) X 2 A 2 , 
the c-quantities from ( |A.1| ) are given in terms of K, h^, and R as 



(A.3) 
(A.4) 



A 



Cll 



C12 



1 + pX 



P + X 

p 1 + pX 
P + X 



CO 



c 2 i 



C22 



1 + z 2 X X 
2 p + X 

l-z 2 XX 
2 p+X 



(A.5) 



C13 = 2Vp^l ^ 

a; 



A 



C 23 = V^T[1+P(1-X 2 )A]-. 



Taking the limit R — > 00 in ( |A.5| ), we get en = C13 = c 2 i = c 2 2 = and 
C12 = c 2 3 = 1. This limit, understood as Eab - ► 00 [see ( j2.2|) l, corresponds 
to the original model introduced by Wheeler and Widom [1|]. 

Let us remark that x is an °dd function of h^, while A and u are even 
functions of [see (A.3) and (A.4)]. As a consequence, the quantities en, 
C12, and c 2 3 are symmetric under the transformation — > — whereas C13, 
c 2 i, and c 2 2 are antisymmetric. Because m is antisymmetric in and cr 
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is an even function of m (see Section 3), from the mole fractions equations 
( 2.12j) to ( 2.14 ) it follows that Xab is symmetric in h+, while Xaa ~ Xbb 



is antisymmetric. These properties determine the symmetry of the phase 
diagrams (discussed in Section 4) under reflections around the (Xaa = X BB )- 
axis. Moreover, at = all Ci 3 = c 2 i = c 2 2 = and o~ 3 _ l2 decouples from 



ni 



3-12 



see ( p!5|) and 



Another interesting limit is when K — > K max [see ( |2.20|) 1. In this case 



A — > oo, uj goes to infinity as A 2 , and from ( |A.5| ) we get Cn = 1, c i2 = c ri = 0; 
consequently, u 3 _ 12 given by ( [2.15|) is equal to one. From (|2.12|) and Q2.13| ) it 



then follows that Xab = 0. This property is used in Section 4 to show that 
in the limit K —>■ K max the coexistence and the spinodal curves always reach 
the bottom side of the triangular diagram. 
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Figure 1: Configuration of molecules within the extended Wheeler- Widom 
model on the honeycomb lattice. 
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Figure 2: A portion of the 3-12 lattice, each site being covered by one C3 
graph (triangle) and one C2 graph (segment connecting two triangles). 
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Figure 3: Star-triangle transformation (up) and double-decoration (down) by 
which the Ising model on the 3-12 lattice can be transformed into the Ising 
model on the honeycomb lattice. At constant L±, the double-decoration can 
be also inverted (i.e., from right to left). 
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Figure 4: The Cayley tree with the coordination number q = 3 and the 
number of shells n = 4. The Bethe lattice corresponds to the limit n — > oo. 
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Figure 5: Isothermal phase diagrams of the extended Wheeler- Widom model 
on the honeycomb lattice when the equivalent Ising model is treated in the 
mean-field approximation (t is the reduced temperature of the model). The 
coexistence curves are drawn by solid lines, while the spinodal curves by 
dashed lines; the solid circles mark the corresponding plait points. At t — 
the two curves (coexistence and spinodal) coincide. 
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Figure 6: The same as in Fig. || but in the Bethe- lattice case. At t = the 
spinodal curve does not coincide with the coexistence curve. 
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Figure 7: The phase diagram of the extended Wheeler- Widom model on 
the honeycomb lattice in the limit of two components (Xab = 0) and when 
the equivalent Ising model is treated in the mean-field approximation. The 
coexistence curve is represented by a solid line, the spinodal by a dashed line, 
and the plait point by a solid circle. 
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Figure 8: The same as in Fig. [?], but in the Bethe-lattice case. 
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